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Transforming  the  Navler-Stokes  equations  for  the 
case  of  turbulent  fluid  motion  by  replacing  the  ac¬ 
tual  parameters  of  the  moving  medium  by  a  sum  of  some 
average  value  and  by  a  pulsation  component  leads  to  an 
extremely  compllr.-.  red  system,  the  solution  which  is 
usually  atterapte.  :y  two  methods:  with  the  aid  of  a 
statistical  ensemble  or  by  using  certain  seml-emplrlcal 
relationships.  The  semi-empirical  theory  of  turbu¬ 
lence  allows  some  results  to  be  obtained,  in  partic¬ 
ular,  the  solution  to  the  problem  of  the  propaga¬ 
tion  of  a  free  tiirbulent  Jet  of  incompressible  fluid 
in  a  quiescent  medium.  However,  this  problem  had 
remained  unsolved  imtll  recently  because  of  the  ab¬ 
sence  of  the  appropriate  seml-emplrlcal  relationships 
needed  to  close  the  system  of  equations  which  describe 
the  propagation  of  a  gas  Jet  in  a  moving  medium.  And 
only  after  the  work  of  G.  N.  Abramovich  [1],  who,  by 
developing  L.  Prandtl*  s  idea,  found  a  seml-ernpirical 
equation  relating  the  intensity  of  the  turbulent  ex¬ 
pansion  of  the  Jet  to  the  values  of  the  flow  velocity 
at  the  boiindaries  of  the  mixing  zone,  were  the  prob¬ 
lem  of  the  propagation  of  a  fluid  Jet  In  a  flow  with 
the  same  physical  properties  and  the  pTOblem  of  gas 
Jets  In  moving  and  quiescent  media  [2-4]  solved. 

A  comparison  between  obtained  solutions  and  the 
experimental  data  of  different  authors  [J]  has  shown 
that  seml-emplrlcal  theory  properly  reflects  the 
general  outlines  of  the  actual  process  of  gas  Jet 
propagation.  The  systematic  discrepancy  observed 
between  the  theoretical  and  experimental  insults  when 
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the  effect  of  the  Initial  turbulence  of  the  inter¬ 
mixing  flows,  which  is  not  accounted  for  in  the 
theory.  However  calculations  show  that  the  ini¬ 
tial  degree  of  this  turbulence  should  be  very  high, 
of  the  order  of  10-15^.  Since  for  steady-state  flow 
in  different  channels  the  degree  of  turbulence  does 
not  usually  exceed  1-5^  (these  same  values  charac¬ 
terize  the  turbulence  of  a  Jet  and  concomitant  flow)^ 
it  is  clear  that  the  abovementioned  qualitative 
explanation  is  gravely  in  need  of  experimental 
grounding. 

Calculation  of  the  acquired  flow  rats  of  a  gas  Jet 
being  propagated  in  a  quiescent  or  moving  medium 
has  shown  that  in  the  region  where  the  theoretical 
and  experimental  characteristics  of  the  Jet  coincide 
the  ejection  properties  of  the  Jets  remain  the  same 
as  for  a  corresponding  submerged  Jet  (with  the  same 
Initial  momentum  and  cross-sectional  area  of  the 
exhaust  nozzle)  of  incompressible  fluid. 

This  fact  suggests  the  idea  of  the  universality  of 
the  ejection  properties  of  a  turbulent  Jet  under  what¬ 
ever  external  conditions  its  growth  occurs.  The 
expressed  hypothesis  permits  us  without  particular 
difficulties  to  use  the  well  known  and  extensively 
tested  theory  of  a  submerged  turbulent  Jet  of  in¬ 
compressible  fluid  to  find  the  laws  for  the  propaga¬ 
tion  of  a  Jet  under  more  complicated  conditions  (non- 
Isothermlc  state,  off-design  efflxix,  dynamic  com¬ 
pressibility,  longitudinal  pressure  gradients,  etc.) . 

We  note  that  this  hypothesis  is  based  on  the 
following  considerations';  since  the  basic  character¬ 
istic  of  a  Jet  is  its  momentum;  and  its  chief  property, 
the  ejection  of  substance  from  the  ambient  medium,  it 
is  natural  to  assume  that  when  the  momentum  of  the  Jet 
is  kept  constant,  its  ejection  properties  remain  \xn- 
changed. 


Analytically,  the  hypothesis  concerning  the  universality  of  the 
ejection  properties  of  turbulent  Jets  is  formulated  in  the  following 
way.  If  a  given  Jet  propagating  under  ordinary  conditions  (subscid.pt 
1)  and  a  standard  submerged  isothermic  free  Jet  ( subscript  a)  issue 
from  Identical  nozzles  and  have  equal  momenta  at  the  initial  cross 
section 

(0.1) 
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then  these  Jets  also  possess  the  same  ejection  properties.  In  other 
words,  the  laws  for  the  increase  in  the  mass  of  a  Jet  due  to  the 
substance  drawn  into  it  from  the  ambient  medium  coincide  for  the 
given  and  standard  Jets  when  condition  (O.l)  is  fulfilled; 


dx  dx 


(0.2) 


Here  G  is  the  gas  flow  rate  in  an  arbitrary  cross  section  of 
the  Jet,  and  x  is  the  longitudinal  coordinate  referred  to  the  nozzle 
edge. 


It  Is  easy  to  show  that  when  condition  (O.l)  is  observed  and 
when  the  areas  of  the  outlet  apertures  are  equal  (besides  having  the 
same  geometric  shape)  the  ratio  of  the  flows  of  the  standard  and 
arbitrary  Jet  Is  expressed  by  the  formula 


c-  =  j /  ^  yjo 

V  So 


(0.3) 


where,  by  assumption,  0  _  Is  the  density  of  the  ambient  medium  Into 
which  the  Jet  issues. 

Let  us  use  formula  (0.3)  to  reduce  relationship  (0.2)  to  dimen¬ 
sionless  form; 


«/*•  *  • 


</*• 


(0.4) 


where  To  Is  the  Initial  radius  of  an  axially  symmetrical  Jet  or  the 
Initial  half-width  of  a  plane-parallel  Jet. 

The  hypothesis  proposed  above  may  be  formulated  in  yet  another 
way;  the  intensity  of  the  accxmiulatlon  of  mass  by  the  Jet  (the  Jet's 
ejection  capacltj)  is  proportional  to  Its  Initial  momentum  dG/M  =  Idem. 

If  there  are  twa  .Jets  .  flowing  out  from  geometrically  similar 
nozzles.  It  Is  not  difficult  to  show  that  the  ratio  of  the  reduced 
increments  (dG“  =  dG/Go)  In  these  Jets  may  be  expressed  by  the  formula 
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When  the  Initial  momenta  and  outlet  aperture  areas  are  Identical 
we  have  dG?  =  “\^^dG|  ,  which  Is  analogous  to  relationship  (0.4)  ob¬ 
tained  above. 

Finally,  It  Is  possible  to  show  that  the  proposed  hypothesis  for 
the  case  of  a  submerged  gas  Jet  conforms  In  the  first  approximation 
to  the  assumption  regarding  the  universality  of  the  profile  relative 
to  the  transverse  component  of  velocity:  =  idem. 

This  assumption  together  with  the  fact,  established  by  many 
researchers  (cf.[3]),  of  the  similarity  of  the  profiles  of  the 
longitudinal  velocity  component  =  Idem)  Is  equivalent 'to  a  hypoth¬ 

esis  regarding  the  universality  of  the  kinematic  flow  pattern  in 
gas  Jets. 

Let  us  illustrate  this  by  the  example  of  an  axially  symmetrical 
Jet  for  which  the  flow  increment  in  an  arbitrary  cross  section  of  the 
submerged  Jet  rfC  =  2ngpoo  Is  represented  In  dimensionless  form: 


=  2— 

dx‘  ~  r,  u»  ^ 


As  will  be  shown  below  [cf.  formula  (2.4)],  the  relationship 

^  _ 1 _ 

r,  ■!«  Yf'O 

Is  valid  for  the  case  under  consideration,  where  the  function  ♦  may 
be  replaced  by  the  constant  O  (u„./ const,  for  not  too  high  values  of 
p • .  Then 

.  1  . 

const  — 

rfx*  Y  f>’  “ 

The  approximate  conformity  between  the  stated  hypothesis 
[cf.  formula  (0.4)]  and  the  assumption  regarding  the  kinematic  uni¬ 
versality  of  the  Jets  v/u^  =  idem  and  u/u^^^  =  idem  also  follow  from 
this. 

Analysis  of  known  experimental  and  theoretical  data  of  various 
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authors  shows  that  the  accumulation  of  mass  along  the  submerged  Jet 
of  Incompressible  fluid  Is  satisfactorily  described  by  the  following 
generalized  relationships: 

for  the  main  segment  of  an  axially  symmetrical  Jet 

=  a, (0.5) 

for  the  main  segment  of  a  plane-parallel  Jet 

c-  (a‘)  =  3;  V  x-  (0.6) 

for  the  Initial  segment  of  a  planar  Jet 

C‘' -f  1  (0.7) 

This  formula  may  also  be  used  for  calculating  the  Initial  portion 
of  an  axially  symmetrical  Jet  in  the  first  approximation.* 

In  talking  about  the  formulas  describing  the  ejection  of  sub¬ 
stance  from  the  ambient  medium  by  a  turbulent  Jet,  It  should  be  noted 
that  they  are  noticeably  dissimilar  for  different  authors  although 
they  are  similar  in  structure  and  In  their  general  shape  satisfy 
relationships  (0.5)“(0.7).  For  example,  for  the  main  segment  of  an 
axially  symmetrical  Jet  of  incompressible  fluid  the  coefficient  Oi, 
according  to  the  data  of  several  authors,  varies  from  0.11  to  0.20, 
althovigh  the  most  frequently  encountered  values  of  Uj,  according  to  the 
data  in  some  papers  [3-9],  He  in  the  Interval  0.13-0.16.  The  main 
reason  for  sucl)  discrepancies  as  well  as  for  the  lack  of  conformity 
among  the  formulas  for  determining  the  attenuation  of  the  axial  velocity 
of  the  Jet  lies,  in  our  opinion,  in  the  carelessness  of  the  analysis 
of  the  experimental  data.  Thus  authors  rarely  present  data  on  the 
distribution  of  velocity  and  other  flow  parameters  in  the  initial 
cross  section  of  the  Jet,  l.e.,  initial  turbulence,  but  relate  the 
axial  velocity  usually  either  to  a  maximum,  or  to  an  average  (with 
respect  to  rate  of  flow  or  momentimi)  value  of  the  velocity  in  the 
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initial  cross  section  besiaes  not  stating  which  of  these  methods  Is 
being  used.  As  will  be  clear  from  what  follows,  taking  into  account 
the  nonuniformity  of  the  flow  at  the  outlet  from  the  nozzle  ( charac¬ 
terized  by  the  coefficients  and  n^)  introduces  Important 

corrections  into  the  calculation,  and  it  should  necessarily  be 
carried  out  when  analyzing  experimental  data  or  when  comparing 
experimental  and  theoretical  results. 

Let  us  present  values  of  the  coefficients  in  foimiulas  (O.p)- 
(0.7)  according  to  available  experimental  data 

a,  =0.i;5,  3, --4.3,  a,  =  0.37G,  ?,  =  0,  -  =  0.0362 

We  note  that  with  the  accu-mulation  of  experimental  material  and 
after  careful  analysis  of  already  known  experimental  data  pertaining 
to  isothermic  submerged  jets  under  fixed  conditions  in  the  Initial 
cross  section.  It  may  become  necessary  to  alter  these  values  some¬ 
what  (in  particular,  the  coefficients  In  formulas  (0.5)-(0.7)  nay  prove 
to  be  dependent  upon  the  initial  degree  of  turbulence  In  the  jet) . 
Therefore  to  preserve  the  generality  of  future  calculations  we  shall 
for  the  time  being  use  the  ejection  laws  in  the  form  of  generalized 
relationships  (0.5)-(0-7). 

It  is  not  difficult  to  obtain  the  ejection  law  for  ar.  arbitrary 
Jet  from  formulas  (0.4)-(0.7).  Thus  we  have  for  the  initial  segment 
of  the  Jet 

c»  =  VT=ti-+i  (0,8) 

while  for  the  main  segment  of  an  axially  synmetrlcal  Jet  we  obtain 

^  (x°  +  3,)  ( ®  •  9) 

In  the  derivation  of  these  relationships  the  constants  of  Inte¬ 
gration  were  determined  from  the  condition  that  the  obtained 


relationships  coincide  with  the  ejection  law  for  an  Isctheralc  sub¬ 
merged  Jet  when  p"  =  1.  Generally  speaking,  these  Integration  con¬ 
stants  may  be  functions  of  the  Initial  parameters  of  the  Jet  (e.g.,p*), 
but  here  we  shall  consider  them  to  be  constants. 

As  may  be  noted,  no  quantities  characterizing  the  kinematics  of 
the  ambient  medium  enter  Into  the  right  sides  of  Eq.  (0.8)  and  (0.9)- 
This  means  that.  In  accordance  with  the  adopted  hypothesis,  the  rela¬ 
tive  velocity  of  the  concomitant  flow  has  no  effect  on  the  value  of 
the  acquired  mass  of  the  Jet.  The  fact  Just  noted  Is  rather  curious 
and  may  be  of  Interest  in  problems  having  to  do  with  mixture  and 
combustion  In  a  stream.  In  particular,  it  follows  from  this  that 
mixing  of  the  Jet  with  the  ambient  medium,  evidenced  by  the  Increase 
In  the  mass  of  gas  drawn  Into  the  Jet,  will  also  occur  when  the  ve¬ 
locities  of  the  Jet  and  concomitant  flow  are  equal.  However,  If  the 
Ecmcntuff.  of  the  ccncomitajit  flew  Is  greater  than  that  of  the  Jet, 
the  flow  will  be  determined  by  the  characteristics  of  the,  concomitant 
flow,  since  In  this  case  the  dominant  role  is  played  by  the  ejection 
capacity  (proportional,  according  to  the  adopted  (-.ypothesls,  to  the 
Initial  momentum)  of  the  concomitant  flew,  rather  than  of  the  Jet. 


Pig.  1.  Flow  pattern  in  mixing  zone  of 
Initial  section  of  Jet. 
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Wie  hypothesis  put  forth  above,  expressed  by  Eq.  (0.4),  Is  an 
alternative  to  the  equation  of  turbulent  gas  Jet  propagation  suggested 
by  the  author  together  with  G.  N.  Abramovich  in  an  earlier  monograph 
[3],  and  together  with  the  equations  of  conservation  of  momentum  and 
excess  heat  content  forms  a  closed  system  for  calculating  a  gas 
Jet  In  a  moving  medlmn. 

In  the  Individual  examples  below  we  shall  show  what  results 
are  obtained  In  the  practical  use  of  the  proposed  hypothesis,  and 
also  a  comparison  will  be  made  between  these  results  and  corresponding 
experimental  data. 

1.  Initial  segment  of  a  Fluid  Jet  In  a  Concomitant  Flow. 

Let  us  consider  the  initial  segment  of  an  isothermic  titrbulent  Jet 
being  propagated  In  a  moving  medium. 

As  shown  In  the  literature  [l],  on  the  basis  of  an  assumption 
of  the  universality  of  dimensionless  excess  velocity  profiles 


3 


It  Is  possible  to  obtain  from  the  equations  of  conservation  of 
momentxim,  flow  rate  and  transverse  fluid  equilibrium  for  the  circuit 
KyiyzOM  (Fig.  1)  the  following  dependence  of  the  relative  ordinate 
of  the  Jet  boundary  layer  on  the  parameter  ra  =  Ua/ui: 

^  =  0.4l6  +  0.134m  (1.1) 

In  order  to  sdlve  the  problem  It  now  remains  to  find  the  varia¬ 
tion  In  the  thickness  of  the  free  boundary  layer  b  along. the  length  of 
the  Jet.  For  this  purpose  G.  H.  Abramovich  has  suggested  [l]  a  rela¬ 
tionship  derived  from  L.  Prandtl's  theory  on  the  mechanism  of 
mixing  In  the  Jet.  But,  as  noted  above,  this  relationship  gives 
Incorrect  results  when  m  >  0.35-  Instead  of  this  relationship  we 
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shall  tiy  to'^make  use  of  our  hypothesis  on  the  universality  of  the 
ejection  properties  of  a  turbulent  Jet.  For  this  we  shall  compute 
the  amount  of  fluid  In  an  arbitrary  transverse  cross  section  of  the 
initial  section  of  a  plane-parallel  Jet 

Vi 

G  =  2gPiUi  (t*  —  y,)  -f  2^'  ^  pudy 
For  this  case  p“  =  1  we  obtain  from  this 

G“=  =  1 (0.554-0.45  (1-2) 

Substituting  into  here  the  function  G^fx®)  according  to 
formulas  (0.8)  when  p®  =  1  and  replacing  the  quantity" y j/b  in 
accordance  with  Eq.  (l.l),  we  obtain 

/,<■  - _ _ 1“ 

U.iyi4-0.31Cm  * 


It  follows  from  this  when  7  =  0.0;562 


1  +  2.3C  »» 


or 


(1.5) 

fl.4) 


Fig.  2.  Dependence  of  reduced  thickening 
coefficient  for  the  mixing  zone  of  an  iso¬ 
thermic  Jet  (p*  »  1)  on  the  parameter  m: 
curve  1  Is' drawn  according  to  Eq.  (l.4]r; 
curve  2,  according  to  Eq.  .(l.5)i  points 
3)  pertain  to  O.  V.  Yakovlevskly ’ s  ex¬ 
periments  [2]j  points  4) to  B.  A.  Zhestkov’s 
experiments  [4], 
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In  Fig.  2. function  (1.4)  is  compared  with  the  corresponding 
experimental  data  obtained  by  the  author  [2]  and  by  B.  A.  Zhestkov 
[4],  as  well  ay  with  a  similar  functional  relationship  obtained  by 
0.  N.  Abramovich  [1]: 


As  follows  from  the  comparison  of  the  theoretical  and  experimental 
results,  Eq.  (1.4),  which  has  been  constructed  on  the  basis  of  the 
adopted  hypothesis,  corresponds  to  the  experimental  points  better 
than  Eq.  (l.5)  In  the  range  0.35  <  m  <  1,  while  when  m  <  0.35  both 
theoretical  curves  agree  satisfactorily  with  the  experimental  values. 

2.  A  Turbulent  Jet  of  Heated  Gas  (with  Variable  Thermodynamic 
Properties)  in  a  Concomitant  Flow.  Consider  the  main  segment  of  a 
turbulent  Jet  of  heated  or  cooled  gas  issuing  into  a  quiescent  or  moving 
raedlvim  of  the  same  composition  as  the  substance  of  the  Jet  (Plg.  3)  , 

We  note  that  recently  in  connection  with  the  development  of  a  differ¬ 
ent  kind  of  system  the  interest  in  Jets  of  Intensely  heated 
( 10, 000-20, 000 ®K)  gas  has  Increased  considerably. 


Pig.  3.  Diagram  of  flow  over  the  main 
segment  of  a  Jet  being  propagated  in  a 
concomitant  flow. 


At  these  temperatures  physicochemical  transformations 
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(dissociation,  ionization)  are  observed  in  the  gas,  which  lead  to  a 
significant  change  in  its  thermodynamic  properties,  particularly  in 
the  heat  capacity,  which  must  not  be  considered  constant  under  these 
conditions  anyway,  as  it  usually  is  in  the  calculation  of  non-iso- 
thermic  Jets  [3].  Analysis  of  prepared  data  on  the  thermodynamic 
properties  of  “Parlous  gases,  such  as  air  [lO],  shows  that  from  0.001 
to  1000  atm.(abs.)  and  from  1000  to  12,000“K  (this  limit  may  apparently 
be  raised  to  20,000°K)  the  equation  of  state  may  be  written  in  the 
form 

(2.1) 


where  1^  is  the  enthalpy,  and  the  exponent  a  has  the  value  0.83  for 
air.  According  to  preliminary  calculations  this  exponent  is  of  the 
same  order  of  magnitude  for  other  gases:  thus  for  hydrogen  a  *  0.92. 
In  the  future  we  will  use  the  equation  of  state  written  in  the  form 
(2.1)  when  P  =  Po,  assuming  that  free  turbulent  Jets  are  Isobaric. 

Experimental  data  on  heated  air  Jets  [3]  allows  us  to  assume 
that  the  relative  excess  velocity  profiles  in  them.  Just  as  in  iso¬ 
thermic  Jets,  are  universal  and  satisfactorily  described  by  the 
following  equation: 


(2.2) 


In  addition,  we  assume  that  the  relative  excess  enthalpy  profiles 
are  universal,  since  at  gas  temperatures  no  higher  than  1000®K  they 
conform,  on  account  of  the  constancy  of  the  specific  heat,  to  the 
relative  excess  temperature  profiles  of  the  Jet,  which,  as  we  know 
[3]  are  tmiversal.  We  will  consider  the  enthalpy  profiles  to  be 
described  by  the  equation 


In  Ollier  to  determine  the  kinematic  and  thermodynamic  parameters 
at  an  arbitrary  point  in  the  Jet  of  heated  gas  taking  into  account 
relationships  (2.2)  and  (2.3)  It  remains  to  find  the  variation  of  the 
axial  velocity  and  enthalpy,  as  well  as  of  the  thickness  of  the  Jet  as 
a  function  of  the  distance  to  its  Initial  cross  section.  For  this 
purpose  we  have  available  a  system  of  three  equations. 

The  first  two  equations  express  the  laws  of  conservation  of 
momentum  and  excess  heat  content  In  the  Jet  and  after  simple  trans- 
formations  (cf.,  for  example,  Abramovich  [3])  may  be  reduced  to  the 
following  form: 
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(2.4) 

(2.5) 


(2.6) 


the  subscript  0  pertains  to  the  initial  cross  section  of  the  Jet; 

00 ,  to  the  ambient  medium;  and  m,  to  the  Jet  axis. 

Let  us  point  out  that  the  equations  written  above  are  of  a  general 
nature  and  only  the  assumption  of  the  universality  of  the  profiles  of 
Au*  and  Al*  was  used  In  their  derivation.  In  addition  these  relation¬ 
ships  are  suitable  for  both  plane-parallel  and  axially  symmetrical 
Jets;  for  these  cases 

df 


_  dy  ^ 
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respectively. 

The  parameters  ^2u‘  ^’^^roduced  above  characterize 

the  nonuniformity  of  the  velocity  and  enthalpy  (or  temperature) 
profile  In  the  Initial  cross  section  of  the  Jet;  In  a  Jet  with 
absolutely  uniform  velocity  and  enthalpy  they  have  the  sairae  value, 
equal  to  unity. 

It  should  be  emphasized  that  the  difference  in  ncnuniformlty 
parameters  from  unity  Is  often  neglected  In  calculations  of  turbulent 
jetsj  as  this  gives  results  which  differ  significantly  from 

actuality.  This  Is  especially  true  for  the  calculation 

of  a  Jet  In  a  concomitant  flow,  since,  as  Is  apparent  from  an  analysis 
of  Eq.  (2.4),  the  effect  of  the  coefficients  n^^  and  n^  (when  they 
are  different  from  unity)  increases  rapidly  as  the  parameter  ra^ 
Increases  from  Its  zero  value. 

As  an  Illustration  we  shall  point  out  that  for  a  steady-state 
turbulent  velocity  profile  In  the  Initial  cross  section  of  the  Jet 
(l.e.,  when  the  thickness  of  the  boundary  layer  Is  equal  to  the  radius 
of  the  outlet  aperture)  satisfying  a  power  function  with  an  exponent 
of  1/7,  the  parameters  n^^^  =  O.815,  n^  =  O.680. 

We  shall  use  Eq.  (0.9)  as  the  third  relationship  to  close  the 
system  of  equations,  l.e.,  from  now  on  we  shall  consider  only  an 
axially  symmetrical  Jet  as  It  Is  the  most  frequently  encountered  In 
practice.  The  method  of  calculation  for  a  plane  parallel  Jet  Is  analo¬ 
gous  and  based  on  relationships  (2. 4) -(2. 5)  taking  Into  account  the 
ejection  law  for  a  plane-parallel  Jet  described  by  formulas  (0.4)  and 
(0.6). 

It  is  not  difficult  to  show  that  for  the  assumptions  regarding 
the  universality  of  the  profiles  of  Au*  and  Al®  Introduced  above  the 


relative  rate  of  gas  flow  in  an  arbitrary  cross  section  of  the  Jet  Is 
determined  by  the  equation 


^  -r  (1  —  '«o)  ^|Au„°I 


(2.7) 


Here 


G  =  ^  pud;?,  C»  =  \  pj'jlF„,  At=\  ~  ~ 

j  •  •'  P|>»n  ^ 

0  u  u 


(2.8) 


By  solving  Eqs.  (2.4),  (2.5)  >  and  (2.7)  together  we  easily 
obtain  the  dependent  variables  In  which  we  are  Interested 
and  F/Fo  as  functions  of  the  relative  rate  of  flow  G",  the  variation 
of  which  as  a  function  of  the  dimensionless  distance  x“  Is  described 
by  Eq.  (0.9) •  We  have 
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(2.9) 

(2.10) 
(2.11) 


As  will  be  clear  from  what  follows,  the  quantity  varies  slightly 
along  the  length  of  the  Jet,  and  when  the  velocity  and  enthalpy  pro¬ 
files  are  uniform  In  the  Initial  cross  section  (n^^  *  n^^  =  =  l) 

this  quantity  has  a  value  close  to  that  for  an  Isothermic  submerged 
Jet  =  0.755  (in  the  case  of  a  plane-parallel  Isothermic  Jet 
kjj^  =  0.86o)  .  Thus  Eqs.  (2.10)  and  (2.1l)  correspond  to  an  almost 
linear  relationship  between  the  excess  axial  enthalpy  and  excess  axial 
velocity  of  the  Jet. 

To  find  the  functions  AUjj^*(x’)  and  Al^®(x*)  It  Is  required  to 
know  the  values  of  A^  and  Since  the  profiles  of  Au®  and  Al® 

In  transverse  cross  sections  are  known  as  functions  of  while  the 
quantity  p/pom  when  (2.1)  Is  used  as  the  equation  of  state  may  be 
represented  In  the  form 
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_p_  _ _ _ 

|l  r  (A- 


(2.12) 


after  Integration  it  is  possible  to  find  the  coefficients  and 
as  functions  of  the  quantities  po  and  ^i^“.  But,  for  the  fractional 
values  of  ct  given  above,  the  integrals  in  Eqs.  (2.6)  do  not  yield 
elementary  functions  when  relationship  (2.12)  is  substituted  into 
them,  and  therefore  numerical  calculations  were  made  for  their 
determination  in  the  range  of  variation  of  the  quantity  ( Po  - 
from  0.75  to  20.0  when  0.5  <  a  <  1.  An  empirical  analysis  of  the 
results  of  these  calculations  has  shown  that  in  said  range  the  values 
of  A^  and  B^  can  with  an  accuracy  sufficient  for  practice  be  computed 
using  the  following  formulas: 

^  0-257  A* 
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_  j _ o.iy.;».*  ( 2 . 15) 

[1 -4-0.700 (r,-l)  A. 

,,  _ 0-<28/.„*  2?  =  _ - 

ll+0.500(;.,-I)A.„,r  ’  ■  U-i-0.730(;.,-l)Ai„T 

Substitution  of  relationships  (2.15)  into  Eq.  (2.9)  by  taking 
formula  (2.IO)  into  account  and  replacement  of  the  quantity  G® 
according  to  formula  (0.9)  using  the  obvious  equality 


P'  =  Po' 


(2.14) 


allows  us  to  obtain  the  relationship  between  the  excess  axial 

velocity  and  the  distance  i-vio  nozzle  in  final  form 

A“«*  »  1.92  v„ 
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The  variation  in  the  enthalpy  along  the  axis  of  the  Jet 
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l3  determined  with  the  aid  of  Eq.  (2.10)  .  The  coefficient 
contained  in  the  theoretical  equations  may  be  determined  in  the  first 
approximation  using  formula  (2.11),  into  which  are  substituted  the 
quantities  Ai,  Bj,  Ba  compuced  from  relationships  (2.15)  using 
the  approximation 


Au'  =  0.753  — 

As  calculations  show,  the  values 
are  extremely  close  to  the  values  of  k 
mation. 


of  de 
.  fc'und 


(2.17) 

temlned  in  this  way 
in  the  second  approxi- 


The  axial  temperature  of  the  jet  for  a  particular  value  of 
A1  *  is  determined  from  tables  of  therr.odynamic  functions  for  the 

d 

given  gas  ( in  the  case  of  air  -  the  tables  cited  at  the  end  of  the 

article  [lO])  as  a  function  T  =■  T  ( 1  )  ,  where 

m  m  nr 


I’m  =  (‘O  —  I'tc)  -r  ( 00  (2. 18) 

The  thickness  of  the  ^et =  r  - ==  in  an  arbitrary  cross 

section  is  found  from  Eq.  (2.4)  using  the  function  AUj^*(x®),  de¬ 
termined  by  Eq.  (2.15). 

Thus  the  equations  obtained  above  permit  the  variation  in  axial 
velocity  and  temperature  and  also  the  thickness  of  the  Jet  to  be 
determined  as  functions  of  the  dimensionless  coordinate  x*  in  terms 
of  the  parameters  of  the  Jet  and  concomitant  flow  in  the  initial 
cross  section.  The  value  of  the  velocity  and  enthalpy  at  an  arbitrary 
point  in  the  Jet  is  determined  from  their  known  values  at  the  slxIb 
with  the  aid  of  Eqs.  (2.2)  and  (2.5). 

In  order  to  appraise  the  effect  of  the  parameter  po,  which 
characterizes  the  difference  in  the  themodynamlc  properties  of  the 
Jet  and  the  ambient  medium,  on  the  propagation  of  the  Jet,  let  us 


consider  the  simplest  case  of  the  escape  of  a  heated  gas  Into  a 
quiescent  medium  (mo  «•  O) .  In  this  case  from  Eq.  (2.15)  when  Mo  ■  0 
we  obtain  the  following  formula  for  the  axial  velocity  of  the  Jet 
after  replacing  the  coefficients  a,  and  by  their  respective 
numerical  values: 


r  I  ^  -l.  12.4v„ 


(2.19) 


where 


i'. 


0.753 


r ' 

"ru  </*•—>)  [  1 


(2.20) 


From  an  analysis  of  the  law  for  the  damping  of  the  axial  velocity 
of  an  Isothermic  Jet  (2.19)  It  follows  that  In  a  Jet  heated  relative 
to  the  ambient  medium  the  axial  velocity  drops  off  considerably  faster 
than  In  an  Isothermic  Jet,  with  the  value  of  u^  for  a  given  value  of 
X*  being  Inversely  proportional  to  the  square  root  of  the  ratio  of 
the  areas  of  the  ambient  medlrim  and  Jet.  As  an  Illustration  the 
calculated  functions  are  presented  In  Pig.  4  for  an  isothermic 

(po  =  1,  or  p"  =  l)  cind  an  intensely  heated  ( Po  =  26,  or  p*  =  15) 

Jet  of  air  escaping  Into  a  quiescent  medium. 

Also  presented  are  experimental  points  characterizing  the  axial 
velocity  In  an  intense3.y  heated  dr  Jet  according  to  data  obtained 
by  V.  Ya.  Bezmenov  and  V.  S.  Borisov.*  As  Is  apparent  from  the  results 
in  Pig.  4,  theoretical  formula  (2. I9)  obtained  on  the  basis  of  the 
hypothesis  regarding  the  universality  of  the  ejection  properties 


♦The  data  of  this  tmpubllshed  work  was  obligingly  made  available 
to  the  author  for  the  preparation  of  the  present  article. 
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of  turbulent  Jets  agrees  satisfactorily  with  the  experimental  data 
right  up  to  very  high  values  of  the  parameter  po. 

The  damping  curve  for  the  excess  axial  temperature  AT®  In  a  non- 

z& 

Isothermic  Jet  is  also  presented  In  Fig.  4,  and  the  corresponding 
experimental  values,  which,  as  can  be  seen,  are  close  to  the  theoreti¬ 
cal,  are  also  shown  there.  This  fact.  In  our  opinion,  attests  Ihe 
negligibly  small  effect  of  radiant  heat  transfer,  which  has  not  been 
taken  Into  account  In  the  present  calculations,  on  the  thermodjTiamlc 
properties  of  Intensely  heated  turbulent  Jets. 

The  theoretical  and  experimental  values  of  the  thickness  of  a  Jet 
of  Intensely  heated  air  are  compared  in  Fig.  and  there  Is  also 
presented  a  curve  characterizing  the  variation  In  the  thickness  of  an 
isothermic  Jet  along  Its  length.  It  is  apparent  that  the  theoretical 
thickness  of  a  heated  Jet  which  conforms  well  to  experimental  data 
markedly  exceeds  the  thickness  of  an  isothermic  Jet  for  the  same  values 
of  X*. 


Pig.  4.  Variation  in  the  axial 
velocity  u^®  along  an  isothermic 

(p°  =  l)  and  a  heated  (p®  =  15) 
air  Jet;  curves  are  calculated; 
small  circles  pertain  to  ex¬ 
periments  performed  by  V.  Ya. 
Bezmenov  and  V.  S.  Borisov 
for  p®  =  15;  variation  In 

axial  temperature  AT®  along 

m  ^ 

the  axis  of  a  heated  (p®  *  15) 
air  Jet;  curve  Is  calculated; 
points  pertain  to  experiments 
at  p*  =  15. 


Fig.  5.  Increase  In  the  thick¬ 
ness  of  an  isothermic  .(P*  -  1) 
and  a  heated  (p“  -  15)  air 
Jet  along  the  x-axls;  curves 
are  calculated,  points  |>eiv 
tain  to  experiments  at  p®  -  15. 
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Fig.  6.  Variation  In  the  velocity  heat 
Ap®  =  along  the  axis  of  heated 

(p®  =  15)  and  Isothermic  ( p'  =  l)  air  jets; 
curves  are  calculated,  points  pertain  to 
experiments  at  p®  =  I5. 

theoretical  and  experimental  data  presented  In  Fig.  5  permits 
us  to  draw  the  conclusion  that,  in  spite  of  the  difference  in  thick¬ 
nesses  of  isothermic  and  heated  submerged  jets,  the  angle  between 
the  boundary  of  the  main  section  of  the  jet  and  the  x-axls,  apparently 
does  not  depend  on  the  parameter  p®,  and  its  tangent  has  a  constant 
value  equal  to  0.22. 

Finally,  In  Fig.  6  there  Is  shown  a  theoretical  curve  character¬ 
izing  the  attenuation  of  the  velocity  head  along  the  axis  of  a  heated 
(p"  »  15)  air  jet,  and  experimental  points  obtained  by  V.  Ya.  Bezmenov 
and  V.  S.  Borisov  are  plotted  there  also.  A  theoretical  curve  for 
the  axial  velocity  head  in  an  isothermic  air  Jet  Is  shown  for  compari¬ 
son  in  the  same  flgtire.  As  we  see,  the  coincidence  between  the 
theoretical  and  experimental  results  is  entirely  satisfactory  for  the 
case  under  consideration  (p®  «  15). 

A  comparison  of  data  concerning  the  velocity  head  at  the  Jet 
axis  for  p*  »  1  and  p*  «  15  does  not  speak  well  for  L.  A.  Vulls' 
hypothesis  [ll]  regarding  the  universality  of  velocity  head  distribu¬ 
tion  in  gas  jets. 

Theoretical  relationship  (2.19)  satisfies  the  experimental  data 


well  also  for  small  values  of  p®.  It  must  be  noted  however  than  when 
p*  ^  3  (l.e.,  when  the  temperature  of  the  Jet  does  not  exceed  1000 •k) 

the  eqfuatlon  of  state  must  be  used  In  Its  usual  form 

A  -  ■C.  li. 

T 

and  the  specific  heat  Cp  must  be  considered  constant.  Then  in  all 
relationships  obtained  above  It  is  necessary  to  set  a  ■  l  amd, 
consequently,  po  -  p*-  Calculations  performed  by  this  method  for  a 
warned  (p*  -  2),  submerged  (iro  =  O)  Jet  are  compared  with  the 
corresponding  experimental  data  [7]  in  Fig.  7.  As  we  see,  a  quite 
distinct  coincidence  between  the  experlm.ental  and  theoretical  results 
is  observed  in  this  case. 

A  comparison  of  curves  for  the  attenuation  of  the  excess  axial 
velocity  of  a  non-lsothermlc  Jet  (p®  *=  2.04)  propagating  in  a  con¬ 
comitant  How  (0<mo<0.5)  with  experlmer.tal  data  [12]  obtained  for  analogous 
conditions  Is  presented  in  Fig.  8,  where  is  taken  ei^  lal  to 

0.95  for  all  values  of  mo.  It  follows  from  an  analysis  of  this  figure 
that  the  hypothesis  regarding  the  universality  of  the  ejection 
properties  of  turbulent  Jets  Is  Justified  in  the  general  case  of  the 
propagation  of  a  non-lsothermlc  Jet  in  a  moving  medium. 


Pig.  7.  A  comparison  of  theoretical  and 
experimental  data  on  the  axial  velocity  of 
a  warmed  up  (p®  ■  2)  and  an  isothemlc 
(p®  »  l)  Jet  in  a  concomitant  flow 
(nj^=  0.888  and  n^^  -  0.8o6)  j  curves  are 

calculated;  points  pertain  to  experiments 
conducted  by  Yu.  V.  Ivanov  and  Kh.  N.  Suy 
[7]  for  the  values;  l)p*  -  2,  mo  ■  0; 

2)p*  -  1,  mo  =  0;  3)p®  »  nvj  -  O.65. 
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Pig.  8.  Varlatlor.  of  the  excess  axial 

velocity  Au®  in  a  non-lscthermic 
ro 

(p®  =  2.04)  Jet  propagating  in  a  con¬ 
comitant  flow  ( ^•95);  curves 

are  calculated,  points  pertain  to  ex¬ 
periments  conducted  by  0.  Pabst  [12]. 


3.  A  Jet  of  Inccmpresslble  Fluid  In  a  Ccnccmitant  Flow.  In 
order  to  more  clearly  Isolate  the  role  of  the  parameter  mo  in  the 
propagation  of  turbulent  Jets,  let  us  consider  the  special  case  of  an 
Isothermic  (p°  »  l)  Jet  issuing  into  a  fluid  flow  parallel  to  it. 

Under  these  conditions  (p®  =  l)  using  formulas  (2.2)  and  (2.6)  we 
obtain 

Ao  =  1,  Ai  =  0.257,  Az  =  0.134  (3.1) 

and  the  final  equation  determining  the  law  for  the  attenuation  of  the 
axial  velocity  assumes  the  form 

(3.2) 

”  3  SJ-'.;.,  -  .ii,  ,  ^  ThTTi 

The  thickness  of  the  Jet  (r®  =  r/ro  -Vp/Fo)  for  a  given  value 
of  X®  is  found  from  equation  (2.4),  the  values  of  the  coefficients 
entering  into  it  are  determined  by  equalities  (3.I). 

In  order  to  Illustrate  the  characteristics  of  the  spreading  of 
an  isothermic  fluid  Jet  in  a  concomitant  flow  (when  »  n^-  l) 

ctirves  of  the  attenuation  of  the  excess  axlail  velocity  Au^  and  of  the 
increase  in  the  thickness  of  the  Jet  r®  when  0  ^  m©  ^  1.  As  also 
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rollowa  from  G.  N.  Abramovich ' s  theory  [3],  the  presence  of  concomitant 
flow  leads  to  an  Increase  In  the  long  range  nature  of  the  Jet,  but  as 
the  parameter  mo  approaches  unity  the  displacement  of  the  axial  ve¬ 
locity  attenuation  curve  becomes  weaker  and  weaker  until  the  curve 
finally  reaches  Its  limiting  position  when  mo  »  1  .  Moreover,  the 
thickness  of  the  Jet  increases  along  its  length  in  this  case.  Just 
as  when  mo  /  1,  whereas  according  to  a  theory  put  forth  in  the 
literature  [3]  the  Jet  does  not  mix  with  the  ambient  medium  when 
mo  ■  1. 

A  comparison  of  the  theoretical  values  of  the  social  velocity  of 
an  isothemlc  jet  in  a  concomitant  flow  ( 0  <  mo  <  l)  with  the 
corresponding  experimental  data  taken  from  the  literature  [7],  is  given 
in  Pig.  7.  Moreover,  we  have  taken  Into  account  the  initial  non- 
uniformity  of  the  Jet,  which  according  to  Yu.  V.  Ivamov's  supplementary 
information,  is  in  keeping  with  n^^  =  0.888  and  n^  *  0.8o6.  As  we 
see,  the  hypothesis  regarding  the  universality  of  the  ejection  prop¬ 
erties  of  turbulent  Jets  is  also  correct  for  the  spreading  of  an 
Isothermic  Jet  in  a  concomitant  flow  the  velocity  of  >4ilch  does  not  ex¬ 
ceed  the  initial  velocity  of  the  Jet. 


Pig.  9»  The  effect  of  the  parataeter 
fflo  on  the  attenuation  of  the  axial 
velocity  of  a  Jet  of  incompressible  fluid 
In  a  concomitant  flow  (n  ■  ”2u"  * 
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Fig.  10.  The  effect  cf  the  parameter 
mo  on  the  thiclcness  of  a  Jet  of  In¬ 
compressible  fluid  In  a  concomitant 
flow  (n^^  =  n^^  =  1)  . 


4 .  A  Supersonic  Jet  in  the  Oesign  and  Cff-I>eslgr;  Efflux, 
ReRlmes.  Let  us  stop  to  consider  the  case  of  a  submerged  supersonic 
Jet  (mo  -  O) .  The  propagation  cf  a  supersonic  Jet  In  a  concomitant 
flow  can  be  analyzed  In  a  similar  way.  The  mcr.entu.m  equation  may  be 
presented  in  the  form  [J] 


"mM:-  -y  r,,,n  -  — -  0  T— , 


(4.1) 


Here  P*  and  P*  are  the  actual  and  design  pressures  In  the 

o 

receiver;  Xo  Is  the  velocity  ratio;  a*  Is  the  critical  sound  velocity; 
and  Cp  and  are  the  specific  heats.  The  cpuantlty  Ag  depends  upon 
the  value  of  aXo*  and  the  axial  velocity  u^  ;  this  functional 
relationship  was  presented  In  the  literature  [3]. 

The  relative  flow  of  the  gas  In  an  arbitrary  cross  section  of 


the  Jet  Is  expressed  by  the  formula 

f, 


(4.2) 


Here  the  coefficient  Ai  Is  a  known  [3]  function  of  aXg  and 
By  combining  relationships  (4.1)  and  (4.2)  we  obtain: 


u  «=  -Ji _ L. 


V  n  J  \  _ 


(4.3) 


t 

If  we  substitute  int^  heru  the  function  G®  -  G”(x*,o“)  according 
to  formula  (0.9)j  where 


(4.4) 


and  Introduce  notation 


At 


(4.5) 


formula  (4.5) 


assumes  the  form 


i.Z)' 


(4.6) 


In  the  first  approximation  the  function  aX©^)  when 

0  ^  aXo®  ^  0.7  and  Uj^^*  ^  1  may  with  an  accuracy  acceptable  for 
practice  be  considered  as  a  constant  S(u^°,  aXo®)  «“  1. 

The  theoretical  curves  for  the  attenuation  of  the  axial  velocity 
1  supersonic  air  Jets  when  aX©-  0.511(Mo  =  1.5)  Is  presented  In 
Pig.  11  for  several  values  of  the  cff-deslgn  parameter  n#  while  In 
Pig.  12  theoretical  curves  are  given  for  design  efflux  with  supersonic 
velocity  when  aX©*  =  0.644  (m©  =  5.0) 


Fig.  11.  Theoretical  and  experimental  data  on  the  velocity 
variation  at  the  axis  of  supersonic  air  Jets  (m©  -  1.5)  for 
design  and  off-deslgn  efflux  regimes;  curves  are  theoretical, 
points  pertain  to  the  experiments  conducted  by  B.  A.  Zhestkov 
et  al.  L3]. 

Pig.  12.  Theoretical  and  experimental  data  on  the  velocity 
variation  at  the  axis  of  supersonic  (m©  ■  5.0,  p°  ■  0.356) 
air  Jets:  curves  are  theoretical;  points  pertain  to 
experiments  conducted  by  B.  X  Zhestkov  et  al.  t3]. 


Also  presented  here  are  the  corresponding  experimental  points  borrowed 
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from  the  papers  published  by  B,  A.  Zhestkov  et.  al.  (cf.[3j)  .  A 
curve  of  the  attenuation  of  axial  velocity  for  an  Isothermic  air  Jet 
with  small  Initial  velocity  (Mo  "•  0,  p®  =  l)  Is  presented  In  Fig.  12 
for  the  sake  of  comparison.  As  we  see  the  supersonic  Jet  Is  of 
considerably  longer  rar.ge,  l.e.,  for  a  given  value  of  x*  Its  axial 
velocity  Is  higher  than  that  of  a  Jet  with  small  Initial  velocity. 
Analysis  of  this  material  Indicates  that  the  calculations  made  on 
the  basis  of  the  hypothesis  about  the  universality  of  the  ejection 
properties  of  turbulent  Jets  also  agree  satisfactorily  with  known 
experimental  data  for  supersonic  (including  off-deslgn)  air  Jets 
right  up  tc  a  gas  velocity  corresponding  to  a  Mach  number  Mo  ”  3- 

5-  Gas  Liquid  Jet.  In  conclusion  we  shall  consider  the  efflux 
of  a  gas  Jet  into  a  medium  filled  with  atomized  liquid.  The  theory 
of  such  a  Jet,  developed  on  the  basis  of  the  author's  formula 

b  I, 

I  [  pudy  (  5 .  l) 

0*  * 

was  first  published  in  Abramovich  [3]* 

We  shall  try  to  apply  our  hy^pothesls  in  this  case  as  well. 

As  has  been  shown  in  the  aforementioned  article  [3],  the  momentum 
equation  leads  to  the  following  relationship  between  the  axial 
velocity  and  the  area  of  a  lateral  cross  section  of  the  gas-liqijld 
Jet; 

(5.2) 

where  N  ■  3.85  for  an  axially  ^mmetrlcal  Jet.  From  the  equation  for  the 
conservation  of  the  initial  amovmt  of  gas  we  obtain  the  following 
relationship  between  the  gas  concentration  at  the  axis  of  the  Jet  and 
the  axial  velocity: 


In  addition.  It  Is  not  difficult  to  show  that  the  flow  rate  of 
the  gas-liquid  mixture  in  an  arbitrary  cross  section  may  be  expressed 
by  the  formula 


C’ 


0.428 

•  A 


(5.4) 


Fig.  13.  A  comparison  of  the  theoretical 
curves  for  the  attenuation  of  the  axial 
velocity  head in  an  air  Jet 
issuing  into  water  (p°  650)  ,  and  In  a 

Jet  with  constant  density  (p“  =  l)  with 
experimental  data  for  a  gas-liquid  Jet 

[13]. 


If  we  substitute  into  here  G®  as  a  function  of  x®  according  to 
formula  (0.9),  relationships  (5.2)  and  (5*4)  form  a  closed  system  of 
equations,  which  permits  us  to  determine  the  velocity  and  concentra¬ 
tion  of  the  mixture  at  the  axis  of  the  Jet,  and  also  the  Jet's 
thickness. 

In  some  cases  the  dynamic  properties  of  a  gas-liquid  Jet,  deter¬ 
mined  by  its  velocity  head,  are  of  interest.  As  has  been  shown  in 
the  literature  [3],  the  velocity  head  of  the  mixture  at  the  axis  of 
the  main  of  the  Jet,  relative  to  the  initial  velocity  head  of 

an  air  Jet  may  be  expressed  by  the  following  approximate  relationship; 


(5.5) 


Substituting  here  u^®  as  a  function  of  x®,  determined  by  relation¬ 
ships  (5.2)  and  (5.4),  we  finally  obtain  for  an  efflux  of  air  into 
water  with  a  high  subsonic  velocity  (p®'*'650) 


-?6- 


A  comparison  of  the  theoretical  curve  for  the  attenuation  of  the 
velocity  head  along  the  axis  of  the  Jet  according  to  formula  (5.6) 
with  experimental  data  borrowed  from  the  literature  [13]  shows 
(Plg.  13)  acceptable  agreement  between  the  method  of  calculation  de¬ 
veloped  here  and  experiment,  especially  If  the  considerable  spread 
of  the  experimental  points  and  the  approximate  nature  of  relationship 
(5-5)  are  taken  Into  account.  The  dotted  line  In  Fig.  13  Joins 
the  curves  for  the  attenuation  of  the  velocity  head  in  the  initial 
(x*  » 2)  and  main  (x®  >  5-3)  segments  of  the  Jet  and  Is  drawn  Ideally. 

A  comparison  of  data  on  the  variation  in  the  velocity  head  along 
the  axis  of  a  gas-liquid  Jet  with  the  corresponding  curve  for  a 
submerged  turbulent  Jet  with  constant  density  (p®  =  l) ,  Just  as  In 
the  case  of  Jets  of  Intensely  heated  air,  shows  L.  A.  Vulls' 
assumption  [ll]  on  the  universality  of  velocity  head  fields  (flow 
of  momentum)  In  turbulent  Jets  to  be  Invalid. 

It  must  be  emphasized  that  In  the  case  of  an  Isothermic  Jet, 
propagating  In  a  concomitant  flow,  the  results  of  a  calculation  based 
on  the  hypothesis  regarding  the  universality  of  ejection  properties 
when  roo  <  0.35  Is  very  close  to  theoretical  data  obtained  on  the  basis 
of  Q.  N.  Abramovich's  theory  [3].  Moreover,  both  the  one  and  the 
other  set  of  results  agrees  well  with  experimental  material.  However, 
at  the  same  time  that  0.  N.  Abramovich's  theory,  for  the  reasons 
pointed  out  above,  begins  to  differ  noticeably  from  the  experimental 
data  as  the  value  of  nvs  approaches  unity,  the  hypothesis  regarding 
the  universality  of  the  Jet's  ejection  properties  entirely  agrees 
with  experiment  even  in  this  range  of  variation  of  the  parameter  Bq. 


We  add  that  the  results  of  calculation  according  to  the  proposed 
method  when  mo  <  0.35  agrees  satisfactorily  with  the  theory  of  gas 
Jets,  developed  by  the  author  and  published  in  monograph  form  [3l- 
This  Is  true  both  for  Jets  issuing  from  the  nozzle  at  high  velocity 
and  for  gas-liquid  Jets.  However,  in  the  tneory  Just  mentioned  non- 
Isothermlc  gas  Jets  are  considered  under  the  assumption  that  the 
specific  heat  Cp  is  constant,  which  is  not  true  when  the  temperature 
of  the  gas  is  considerable  (c®  >  3) •  If  the  dependence  of  Cp  on 
temperature  Is  taken  into  account  in  the  theoretical  formulas,  as 
was  done  in  the  present  work,  the  above-noted  cori'ormity  between  the 
calculations  using  the  new  method  and  using  the  gar,  Jet  theory  [3] 
in  which  the  widening  of  the  Jet  is  determined  from  the  relationship 

Will  probably  also  be  observed  in  the  case  of  strongly  non-lsothermlc 
gas  Jets. 

If  we  take  Into  account  the  fact  that  the  hypothesis  regarding 
the  universality  of  ejection  properties  agrees  satisfactorily  with 
the  results  of  experimental  investigations  of  Jets  with  a 
essentially  variable  density  (plasma,  gas-llquld,  and  supersonic 
Jets)  then  the  method  of  calculation  based  on  this  hypothesis  may 
be  recognized  as  being  sufficiently  reliable. 

Analogous  problems  for  plane-parallel  gas  Jets  may  also  be 
solved  by  the  method  proposed  here.  In  addition  the  method  of  calcu¬ 
lation  based  on  the  hypothesis  regarding  the  universality  of  the 
ejection  properties  of  turbulent  gas  Jets  will  appajoently  also  prove 
applicable  to  the  case  of  Jets  propagating  in  a  space  bounded  by 
solid  walls  (e.g.,  in  the  mixing  chamber  of  a  Jet  pump)  in  analyzing 
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gas  Jets  growing  in  a  drifting  side  flow  as  well  as  In  examining  the 
problems  of  turbulent  combustion  In  a  flow. 

In  conclusion  the  author  takes  the  opportunity  to  thank  0.  N. 
Abramovich  for  his  extremely  useful  advice  in  the  completion  of  this 
work,  and  also  to  express  his  gratitude  to  Ye.  Ya.  Flrsova  who  parti¬ 
cipated  in  the  calculations. 


Received  July  21,  i960. 


RSFERENCES 


1.  Q ■  N.  Abramovich.  0  turbulentnom  smeshenil  na  granitse  dvukh 
ploskoparal'lel'nikh  polokov  zhidkostl  ( prl  sputnem  1  vstrechnom 
dvlzhenli) ,  Sb.  "Teoretlcheskaya  gldromekhanlka" ,  No,  19,  pod  red. 

L.  I.  Sedova.  Oboronglz,  1956. 

2.  0.  V.  YakovlevsVl y ,  X  veprosy  o  tolshchins  cony  tu.rt'iT  entnogo 
peremeshlvaniya  na  granitse  dvukh  petokov  gaza  raznoy  skorosti  1 
plotncstl,  Izv.  AN  SSSR,  OTN,  No.  10,  1958. 

3.  G.  N.  Abramovich.  Teoriya  turbulentnykh  stray',  Fizratgiz, 

i960. 

4.  G.  N.  Abramovich.  Turbulentnaya  struya  v  dvizhushcheysya 
srede,  Izv.  A>4  SSSR,  OTN,  No.  6,  1957. 

5.  G.  N.  Abramovich.  Turbulentnyye  svebodny-.-e  strui  zhldkoatey 
1  gazov,  Gosenex'golzdaFJ  1948. 

6.  V.  N.  Taliyev.  Aerodlnamlka  ventllyatsll,  Stroylzdat,  1954. 

7.  Yu.  V.  Ivyiov  i  Kh.  N.  Suy.  Issledovanlye  razvltlya  strui 
V  sputnom  potoke,  Izv.  AN  EstSSR,  seriya  tekhn.  1  flz.-mat.  n.,  7, 

No.  2,  1958. 

8.  M.  L.  Al^rtson,  J.  B.  Day  and  R.  A.  Jensen.  Diffusion  of 
submerged  Jets.  froc.  of  the  ASCe,  194b,  v.  74,  pp.  1571-I596. 

9.  J.P.  Taylor,  H.  L.  Grlmnett  and  E.  W.  Comings.  Isothermal 
free  Jets  of  air  mixing  with  air.  Chem.  Eng.  Progr. ,  I$51,  v.  74, 

No.  4,  pp.  175-180. 

10 .  A.  S.  P^dyoditelev,  Ye.  V.  Stupochenko,  Ye.  V.  Samuylov. 

I.  P.  Stakhanov,  A.  S.  Pleshanov  and  I.  B.  Rozhdestvenskly.  Tablltsy 
termodlnamlchesKlkh  funktsly  vozdukha  v  Intervale  temperatur  ot 
1,000  do  12,000“K  1  intem/ale  davlenly  ot  0.001  do  lOOO  atm.,  Izd-vo 
AN  SSSR,  1957. 


11.  L.  A.  Vulls.  0  lurbul'.-ntnjkh  f;;azovykh  sti'uyakhj  Izv. 

AN  KazSSR,  serlya  cnergetich. ,  19^^,  vyp.o. 

12.  0.  Pabst.  Die  Ausbreitung  hoisser  Gaastrahlen  in  bewegter 
Luft,  II  Teil.  Deutsche  Luftfahrforschung,  U.  u.  M.  8007,  19^^. 

13.  B.  F.  Gllkman.  0  struye  gaza  v  zhidkostl,  Izv.  AN  SSSR, 
OTN,  Energetika  1  avtomatlka.  No.  2,  2959- 


MIXING  TUR3ULSNT  JETS  CF  DIFFEF^tiT  EENSTTY 

G.  N.  Abranovlch 

A  number  of  papers  [1-h]  have  dealt  with  the  problem  of  the  turbu¬ 
lent  mixing  of  jets  of  different  density.  However,  the  contradictory 
results  of  the  theoretical  papers  of  different  authors  and  the  signif¬ 
icant  discrepancies  In  a  number  of  cases  between  the  theory  and  experi¬ 
ments  force  us  to  return  again  to  this  problem.  Until  recently  the 
following  two  cardinal  questions  in  jet  theory  had  remained  vague: 

1)  Are  the  same  dimensionless  velocity  profiles  retained  in  a 
gas  Jet  as  in  a  Jet  of  incompressible  fluid? 

2)  How  does  one  density  field  or  another  Influence  the  shape  of 
the  Jet  boundary? 

The  material  presented  in  the  author's  monograph  [3]  and  in 
articles  by  Vulis  [4],  Yershin  [5],  and  Gllkman  [6]  show  that  the  first 
qiwstlon  may  be  answered  in  the  affirmative,  l.e.,  it  is  possible  to 
consider  the  velocity  profile  in  different  cross  sections  of  a  Jet  of 
essentially  variable  density  (including  a  combustion  torch  and  a  two- 
phase  gas-liquid  Jet)  to  be  expressed  by  the  same  universal  law  as  in 
the  case  of  an  Incon^jresslble  fluid. 


/' 

The  same  may  be  said  •  •'•nlii,-  u./-:  urai  versallty  of  tlie  profile 
of  the  stagnation  tcmpei'atia-c  and  additive  concentration  for  a  Jet  of 
moderate  velocity  (subsonic  or  low  supersonic)  and  for  a  moderate 
temperature  difference;  .’ic  pi-oblem  of  the  temperature  profile  in  a 
Jet  of  high  supersonic  o:  whe.i  tiic-  dl  f  ft- ''.ri  .in  the  inj>er  a- 

tures  at  the  boundarioa  of  the  mixing  zone  is  great  requli-es  further 
experimental  study. 

The  results  pertaining  t'  ;:ho  second  questli-.’.  pi.c-d  abovr,-  ;u 

less  distinct. 

In  the  author's  paper  [l]  uhe  tuicicness  of  the  mixing- zone  of 
two  Jets  of  different  density  was  evaluated  on  the  assump¬ 
tion  that  the  degree  of  turbulence  in  a  compressible  gas  is  character¬ 
ized  by  the  ratio  of  the  velocity  head  of  the  pulsation  velocity  to 
the  velocity  head  of  the  average  velocity. 

Subsequent  experiments  have  shown  that  the  evaluation  of  the 
Influence  of  compressibility  on  the  thickness  of  the  Jet  obtained  In 
this  case  Is  not  even  reliable  with  respect  to  the  sign  in  a  number 
of  cases  (for  example,  according  to  this  theory  It  luirns  out  that  as 
the  density  increases  in  a  submerged  Jet  the  Jet's  thickness  3-ncreases, 
while  the  opposite  effect  is  observed  in  experiments  [5]) . 

In  0.  V.  Yakovlevskl^/' s  paper  [2]  another  method  was  suggested 
for  taking  Into  account  olic  iiii  luence  of  the  compressibility  of  the 
medium  on  the  thickness  of  the  Jet,  which  method  agrees  satisfactorily 
with  experimental  data  when  the  degree  of  compressibility  for  submerged 
Jets  is  not  very  great,  and  also  In  cases  when  the  velocity  heads  in 
the  mixing  Jets  are  considerably  different.  In  the  authors  monograph 
[3]  0.  V.  Yakovlevskly' s  method  was  extended  for  large  density  gradi¬ 
ents.  ■ 
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However  this  method  of  taking  Into  account  the  conpreaslblllty 
also  proves  not  accurate  enou^  for  very  high  degrees  of  compress- 
Iblllty  (plasma  Jets)  and  altogether  unsuitable  for  the  mixing  of 
Jets  In  which  the  velocity  heads  are  commensurate. 

It  follows  from  theoretical  considerations  that  the  Instability 
of  the  tangential  discontinuity  surface,  leading  to  turbulence,  is 
due  to  the  difference  In  the  velocity  heads  cn  either  side  of  this 
surface,  l.e.,  when  the  velocity  heads  are  the  ^arra  the  discontinuity 
surface  is  most  stable  and,  consequently ,  the  rrlxing  zone  has  minimum 
thickness  In  this  case. 

Experiments  performed  by  L.  A.  Vulls  [4]  and  I.  B.  Palatnlk  [7] 
bear  out  these  considerations  by  shewing  that  regardless  of  the  dif¬ 
ference  In  the  densities  and  velocities  of  the  adjacent  Jets  their 
least  deterioration  Is  observed  when  the  velocity  heads  are  equal, 

O.  V.  Yakovlevskly ' s  method  of  accounting  for  compressibility  [2]  Is 
found  to  be  in  contradiction  with  these  results. 

Proposed  below  is  a  new  method  for  taking  Into  account  the  Influ¬ 
ence  of  the  conpresslblllty  on  the  thickness  of  Jet  mixing  zones,  which 
does  not  contradict  the  experimental  data  for  any  of  the  possible 
mixing  regimes. 

As  in  all  the  above-mentioned  methods,  we  will  consider  the  rate 
of  Increase  In  the  mixing  zone  along  Its  length  to  be  proportional  to 
the  degree  of  flow  turbulence 

(1)  i 

dx  [  u\  1 

I 

I 

Furthermore,  we  may  assume  that  the  momentum  acquired  (or  lost)  j 

by  the  Instantaneous  mass  of  a  Jet  in  connection  with  turbulent  pulsa¬ 
tions  Is  proportional  to  the  difference  In  the  velocity  heads  at  a 

^  ,  ■ 


distance  equal  to  the  mixing  length 


(2) 


Here  pu  Is  the  local  value  of  the  flow  density  of  the  forward 
(along  the  x-axls)  stream;  jv«  j,  the  average  value  of  the  lateral 
pulsation  velocity  (along  the  y-axls) ;  I,  the  mean  value  of  the  mixing 
length;  and  the  thickness  of  the  mixing  zone  at  a  dlstauice  x  from 
the  beginning  of  mixing. 

But  the  lateral  velocity  head  gradient  Is  proportional  to  the 
difference  in  the  velocity  heads  at  the  boundaries  of  the  mixing  zone 
and  inversely  proportional  to  the  zone's  thickness; 

p, 

Hence 


(3) 


Subscripts  1  and  2  pertain  to  the  two  mixing  Jets, 

It  is  necessary  to  give  a  definition  of  the  characteristic  values 
of  the  flow  density  (pu)  and  velocity  (u),  which  must  be  substituted 
Into  the  denominator  of  the  right  side  of  the  last  expression.  We 
shall  assxune  that  some  velocity,  which  is  averaged  for  the  zone 
(characteristic  velocity)  Is  defined  In  the  following  wayi 

{nH) 

while  the  characteristic  values  of  the  density  of  the  flow  and  the 
density  of  the  medium  are  obtained  by  averaging  these  quantities  with 
respect  to  the  thickness  of  the  mixing  tone 

•  • 

Then  If  the  relative  values  of  the  mixing  length  are  universal 
in  different  cross  sections  of  the  mixing  zone  (1/b)  *  idem)  we  have 
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t  =  A(o..ir  5  /  (\  j' 
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The  values  of  the  Integrals  are  calculated  for  given  velocity  and 
density  profiles  [3].  The  four  was  Inserted  Into  the  denominator  of 
the  right  side  so  that  the  coefficient  of  proportionality  (c)  would 
have  the  same  value  as  In  the  author* s  monograph  [3]  for  the  special 
cases  of  the  mixing  of  Jets  of  Incompressible  fluid  and  for  submerged 
gas  Jets.  If  the  density  In  the  Jets  does  not  differ  greatly  (by  less 
than  a  factor  of  two) ,  then  It  Is  possible  to  assume 

Then  ®  « 


<ll>  C  Ojl/I*  —  P,U;* 

i  ^  (Pi‘'i  "f 


(Pi  -r  Ps)  = 


c 


1  — 
(1  + 


(5) 


Here 


•  u* 


_P3. 


In  a  Jet  of  Incompressible  fluid  (p°  =  1)  this  new  expression 
reduces  to  the  expression  previously  [3]  recommended  by  the  author 

db  ^  1  —  m 

dx  ^  1  T-  m 

In  the  case  of  concomltcint  Jets  of  equal  velocity  (m  =  l)  we 
have  from  (3) 

db  _  e  i  —  p* 
dx  ~  2  iT7 

For  a  submerged  Jet  (m  =  O)  when  the  degree  of  compressibility 
Is  not  very  great,  we  obtain  from  (5)  the  same  expression  as  was 
obtained  In  the  literature  [2] 


rT::-7T-63-???/l-2^ 


However,  when  the  difference  In  the  density  of  the  Jet  and  the 
ambient  medium  Is  large  the  new  expression  for  the  thickness  of  the 
mixing  zone  differs  significantly  from  the  old  one  [5],  For  a  sub- 
merged  Jet  (m  =  O)  it  follows  from  W 

b  b 

fih  _  c  r  p  </»/  /  /I*  ?  </*/  \* 

t/x  ”  J  p,  6  /  I  J  Pi  lii  b  ) 

a  o 

whereas  the  previous  expression  had  the  form 

b  b 

c'o  _  <  c  p  Off  /  \  p  ^ 

</x  ~  2  J  pi'  t  /  '  >1  u,  T" 

U  •* 

The  Integ^^als  In  both  of  these  expressions  have  been  computed  for 
any  value  of  p°  In  chapter  7  of  the  author’s  above-mentioned  monograph 
[3].  For  exan?)le,  when  p°  =  10  the  mixing  zone  according  to  new 
expression  (6)  tttrns  out  to  be  almost  twice  as  thick  as  by  former 
expression  (7),  which  agrees  with  the  experimental  data, 

I 

For  equal  velocity  heads  (m®p°  =  l)  there  Is  no  mixing  (db/dx  =  O) 

In  accordeince  with  (5)  .  This  can  be  true  only  when  the  Initial  turbu¬ 
lence  In  the  Jets  Is  small;  If  such  turbulence  is  sufficiently  high, 

I 

then  for  values  of  m*p°  close  to  unity  the  mixing  Is  determined  not  j 

I 

by  the  difference  In  the  velocity  heads  of  the  Jets  but  by  their  ini¬ 
tial  turbulence.  The  proposed  new  method  for  accounting  for  the  effect 
of  the  con5)resslbillty  of  the  medliim  on  the  thickness  of  the  mixing 

zone  requires  careful  experimental  verification  and,  possibly,  fxirther  | 

'  i 

refinement.  However  the  fact  that  this  method  agrees  with  all  avail-  | 
able  experimental  data  In  the  limiting  cases  of  a  submerged  gas  Jet, 
of  mixing  of  Jets  of  Incompressible  fluid,  and  of  mixing  of  Jets  of 
different  density  but  ne£u?ly  equal  velocity  heads  speakS  well  for  It, 
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